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A system consisting of two neutral spin 1/2 particles is analyzed for two magnetic field per- 
turbations: 1) an inhomogeneous magnetic field over all space, and 2) external fields over a half 
space containing only one of the particles. The field is chosen to point from one particle to the 
other, which results in essentially a one-dimensional problem. A number of interesting features are 
revealed for the first case: the singlet, which has zero potential energy in the unperturbed case, 
remains unstable in the perturbing field. The spin zero component of the triplet evolves into a 
bound state with a double well potential, with the possibility of tunneling. Superposition states 
can be constructed which oscillate between entangled and unentangled states. For the second case, 
we show that changes in the magnetic field around one particle affect measurements of the spin of 
the entangled particle not in the magnetic field nonlocally. By using protective measurements, we 
show it is possible in principle to establish a nonlocal interaction using the two particles, provided 
the dipole-dipole potential energy does not vanish and is comparable to the potential energy of the 
particle in the external field. 
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I. INTRODUCTION 



It is of fundamental interest to investigate the prop- 
erties of elementary systems in quantum physics. We 
have considered a system consisting of two neutral spin 
1/2 particles that interact through the spin-spin poten- 
tial. This distinctly quantum system displays many of 
the characteristic features of quantum systems, including 
entanglement, tunneling, bound states, decaying states, 
and spontaneous symmetry breaking. A similar system, 
consisting of harmonically trapped alkali ions that may 
be entangled through the dipole-dipole potential, has 
been proposed for use in quantum computing 0, 0, 0- 
Hence there is practical as well as fundamental interest 
in the model we consider. We wanted to investigate a 
system of two entangled particles that interact through 
a potential that vanishes at infinity. One purpose of the 
model is to allow us to investigate the behavior of one 
portion of an entangled system when we adiabatically 
perturb the other portion, and thereby to investigate the 
coupling between the separate portions of the entangled 
system as a function of the distance between them. As 
mentioned, the properties exhibited may find application 
in quantum computing. In our model, we find a coupling 
between portions of an entangled state that in principle 
allows one to send signals by the modulation of the mag- 
netic field provided the potential energy from the dipole 
potential does not vanish. The maximum separation pos- 
sible is probably of the order of micrometers or less and 
depends on the maximum modulation frequency of the 
signal. In principle, protective measurements jj, |5j can 
be done in one region of the system to determine the el- 
ements of the reduced density matrix, something which 
would not be possible using conventional measurements 



unless an ensemble of identical systems was available 6J . 
The use of protective measurements would permit the 
entanglement to remain. The adiabatic perturbation in 
the protective measurement can be as large as desired, 
so long as the state evolves continuously and the instan- 
taneous energy eigenvalue does not cross that of other 
levels. 



Other interesting features of this spin-spin coupling 
model are apparent when we apply an adiabatic pertur- 
bation that is an inhomogeneous magnetic field over all 
space. We find the initial singlet evolves into an unbound 
state and the triplet develops a double hump potential, 
suggestive of spontaneous symmetry breaking. For the 
latter case, it is possible for one particle to tunnel across 
the barrier to the other side. Also, we are able to study a 
system in which a superposition evolves continuously in 
time, with the wavefunction changing from entangled to 
unentangled and back to entangled. When we eliminate 
a spatial cutoff and allow the dipole-dipole potential to 
become infinite, spontaneous symmetry breaking occurs 
in the degenerate ground state. 

The paper is organized as follows. A model for entan- 
gled states via spin-spin interaction is constructed in Sec. 
II. We study this model for an inhomogeneous magnetic 
field over all space in Sec. Ill, and for a constant and in- 
homogeneous magnetic field over a half space containing 
only one of the particles in Sec. IV. We discuss the possi- 
bility of the protective measurement in our model in Sec. 
V. We summarize our results in Sec. VI. Two appendices 
include another possible model and detail calculations 
outlined in the text. 
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II. THE MODEL 

We will assume that we have a pair (designated 1 and 
2) of identical, uncharged, spin 1/2 particles with coor- 
dinates X\ and x 2 (and corresponding momenta p\ and 
p 2 ). We will apply a magnetic field and determine the 
evolution of the system in the impulsive approximation, 
in which we assume the kinetic energy of the system does 
not change as we apply the magnetic field. We will con- 
sider two cases. Case 1: an inhomogeneous magnetic 
field is present throughout all space, and Case 2: the 
magnetic field is present only in the region to the right 
of the origin. 

For both cases, the Hamiltonian without external fields 
for our system is 



H = p- + ^ + U(\ Xl -x 2 \) 
2m 2m 



(1) 



and the potential energy U for two interacting magnetic 
dipoles fieri and \xcr 2 located at Xi and x 2 respectively 
is 

U = n z : — ^ - -n z S{x 1 -x 2 ) (2) 

Ixi-x^ 3 

where n is a unit vector in the direction (xi — x 2 ). The 
first term is the usual dipole-dipole interaction while the 
second term is the hyperfine interaction term. 

In our paper, we will first analyze the model with the 
approximation that the changes in the kinetic energy, po- 
tential energy, and the relative position of the particles 
are all negligible when we turn on a perturbing magnetic 
field (impulsive approximation) ■ In this approxima- 
tion, the position of the particle does not change signif- 
icantly during the interaction. This approximation is in 
the same spirit as the Born-Oppenheimer approximation 
in which the electronic motion about the nuclei of a di- 
atomic molecule is much more rapid than the vibrational 
motion of the nuclei. Thus it is possible to obtain the 
eigenfunction for the nuclear motion using the energy 
eigenvalue for the electronic motion as the potential. On 
the other hand we assume that during the application of 
the magnetic field the spins evolve adiabatically. This 
approximation is based on the observation that the spin 
precession for the states is much faster than the transla- 
tional motion. With this impulsive approximation, there 
is no change in the potential energy or the relative posi- 
tion of the particles and in this sense the states act as if 
bound. To determine if the states are in fact bound, one 
would have to treat the separation as a dynamical vari- 
able, and include the potential and kinetic energy terms 
in solving Schrodinger's equation. If the potential energy 
eigenvalues obtained with the impulsive approximation 
are positive, then it is very unlikely that the state is in 
fact bound. Only states with negative potential energy 
eigenvalues could be bound. 

Under these assumptions we can choose the coordi- 
nates system so that the particles are separated along the 



z-direction, and z\ and z 2 have opposite signs. There- 
fore we can rotate the coordinate system such that 
(xi — x 2 ) — ► (zi — z 2 )n where n is a unit vector in the z- 
direction. With this coordinate transformations the term 
<Xj • n(i = 1,2) becomes a Zi . Using cr\ ■ <j 2 = 2S 2 — 3 
and neglecting the hyperfine interaction term which is 
relevant only at short distances (see the discussion in 
Appendix A), we can approximate U as 



U(\z\) = ^ 



(2S 2 -3)~3ct Zi <j z 

M 3 



z\ - z 2 . 



(3) 



The classical dipole force will always be in the z-direction, 
and therefore the separation will always be along the z- 
axis, and we have essentially a one dimensional problem. 
For two classical dipoles oriented along the z-axis, the 
interaction would result in an attractive (repulsive) force 
if the dipoles were parallel (antiparallel). 

We will study the effects of applying an adiabatically 
perturbing magnetic field also in the z-direction. With 
this special choice of field, the problem remains a one 
dimensional problem since the magnetic force on each 
particle is also in the z-direction. 

The spin components of the eigenstates of Hq will be 
simultaneous eigenstates of the total spin S and the to- 
tal spin in the z-direction, S z = S Zl + S Z2 . These states 
comprise the usual singlet state with S — 0, S z = 
corresponding to \S), and the triplet spin eigenstates 
with S = 1;S Z = — 1,0, +1 corresponding to \Tg x ) = 
|T_i) , |T) , |7i)). Thus the dipole-dipole potential results 
in entangled statesQ- Because of the indistinguishabil- 
ity of the particles, it is not possible to describe which 
particle is on the left or right, instead quantum mechan- 
ics indicates there is a superposition of both. The spatial 
part of the wavefunctions is chosen so the total wavefunc- 
tion is antisymmetric with respect to the interchange of 
particles 1 and 2. Introducing the symmetrized (anti- 
symmetrized) wave function by 

ip±{zi,z 2 ) = -^(iPr{zi)iPl{z2) ± ^l{zi)^r{z 2 )) (4) 

where ipFi{zi) (i = 1,2) represents the wavefunction for 
particle i on the right side of the origin (z, > 0), and 
^h{zi) on the left side (zj < 0), the singlet state is 

|5) = ^(l+>i|-> a -|-)il+>2)^+(*i,%). (5) 

The spin 1 triplet states for S z = 0, +1,— 1 respectively 
are 

|7 1 > = ^(l+>i|-> 2 + |->il+) 2 )^(^i^2) (6) 



\Ti) = \+) 1 \+) 2 ^-(z 1 ,z 2 ) 



\T- 1 ) = \-) 1 \-) 2 ^{z 1 ,z 2 ). 



(7) 



(8) 
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Using the following relations 

tri-a2\S) = -3\S) (9) 

tri'* 2 \Ti) = \Ti) for alii (10) 

\S) = - \S) (11) 

a Zl a Z2 \T) = -\T) (12) 

^i^s \T±i) = \T±i) , (13) 

the spin-spin interaction U can be represented with the 
basis states {|T_i) , \Ti) , |T) , IS 1 )} as follows: 



the requirement that the divergence of the magnetic field 
vanish 8J. To meet this requirement for an inhomoge- 
neous field in the ^-direction, we need to also have a 
large constant magnetic field in the z-direction. In or- 
der to determine the evolution of the initial singlet or 
triplet state when the magnetic field is applied, we find 
it convenient to decompose the interaction term as 



Hj = -/i(E_+£ + ), 



(18) 



where we define 



U = f{r) 



where we define 



tt \ V 
f( r ) = "a . 



-2 
0-200 
4 




\Zi - z 2 \ 



(14) 



(15) 



The unperturbed potential energies corresponding to 
the eigenstates |T_ X ) , \T X ) , \T) , \S) are -2/, -2/, 4/, 
respectively. The states |T_i) , |Ti) could be considered 
as analogous to classical systems in which the two mag- 
netic moments are parallel to each other, leading to at- 
tractive forces between the particles, whereas the \T) , \S) 
states could be considered analogous to classical systems 
in which the magnetic moments are antiparallel resulting 
in repulsive forces. It is interesting that the singlet state 
has zero energy for the spin-spin interaction and therefore 
is not expected to be a bound state. There is no classical 
analog to this unique quantum mechanical result of zero 
potential energy for the singlet which follows from the 
properties of quantized angular momentum. The triplet 
states characterized by a negative energy are the only 
states that might be bound if the kinetic energy were 
included. In any event, the separation will not change 
significantly as we apply the magnetic field perturbation. 

At time t — we assume we turn on the interaction 
Hamiltonian Hj: 



Hi = -fj,(a zl B(zi,t) + a Z2 B(z 2 ,t)). 



(16) 



We assume that we turn on the B field slowly (adiabat- 
ically) in the z-direction, so that the spin system can 
adjust to the new field and therefore remain in an eigen- 
state of the instantaneous Hamiltonian (adiabatic theo- 
rem) [7j . The total reduced Hamiltonian H R (neglecting 
kinetic energy terms) is now 



H 1 



2 <Tl • <T 2 — 3(T Zl <7 Z2 



^_ _ — _ fJ/ ^ (7zi B(z 1 ,t) + a Z2 B(z 2 , *)) 

(17) 

where B(z\,t) is the value of the magnetic field, which 
always points along the z-direction, at {z\,t). As noted 
previously, this choice of field is done for simplicity and 
is not the most general field that can be applied. We 
also note that in general this field is not consistent with 



£_ = -(a Zl -a Z2 )(B(z 1 ,t)-B(z 2 ,t)) (19) 

£+ = \{a Zl +a Z2 ){B{ Zl ,t) + B{z 2 ,t)). (20) 

There are two cases of magnetic fields that we will con- 
sider. In the first case, an inhomogeneous B field is pro- 
portional to z everywhere. In the second case the ex- 
ternal magnetic field is present on the right side of the 
origin only (zi, z 2 > 0). 



III. CASE 1: INHOMOGENEOUS MAGNETIC 
FIELD IN THE ^-DIRECTION PRESENT IN ALL 
SPACE 

The time independent magnetic field is defined by 



B{ Zi ) = B Q + bZi, i= 1,2 



(21) 



where Bq and b are constants satisfying the condition 
mentioned earlier. For this choice of the magnetic field 
we first apply the large constant field impulsively so that 
there will be no transitions among the states. Then we 
turn on the inhomogeneous part adiabatically [3 . In this 
case we can show that the total Hamiltonian Ht = Hq + 
Hi is separable into two commuting terms that depend, 
respectively, on the center of mass coordinate Z and the 
relative position coordinate z. We define the center of 
mass coordinate Z = \{z\ + z 2 ), with conjugate center of 
mass momentum P = p Zl + p Z2 and the relative position 
coordinate z = z± — z 2l with conjugate momentum p = 
Wz x ~Pz 2 )- Since [p, Z] = [P, z] = and [P, Z\ = [p, z] = 
—i, we obtain the Hamiltonian (the reduced mass is m/2 
for two identical particles) 



p2 2 



(22) 



Using the decomposition ljl8|l the total Hamiltonian is 
rewritten 



H T = H 1 (P,Z) + H 2 (p,z) 



(23) 



where the term depending on the center of mass coordi- 
nate Z is 

p2 

Hi = ^ - MKi + ff« ) (Bo + bZ) (24) 
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and the term depending on the relative coordinate z is 

H 2 = ^- + U(\z\)-^(a Zl -a Z2 )z. (25) 

In the subspace spanned by IS*) and \T) we can show 
that [cti ■ cr 2 , a zi + a Z2 ] = 0. The commutators with only 
terms in a Zl and a Z2 vanish, so [Hi,H 2 ] = 0, and the 
total energy is the sum of the energy eigenvalues for Hi 
and H 2 . 

First we find that the total reduced Hamiltonian H R 
is expressed in the basis {|T_i) , |Ti) , |T) , |5» 



If the magnetic field, which is proportional to b, van- 
ishes, then sinu> — ► and \—a) — ► |S), |+a) — ► |T), as 
expected. For \—a), the effective potential E R (z) is neg- 
ative and always concave down so no stable singlet states 
are expected (Figures 1 and 2). For \+a) the effective po- 
tential E R (z) is a positive double hump potential which 
suggests stable states as shown in Figures 3 and 4. The 
peak near z\ w z 2 ~ is due to the rapid increase in the 
dipole-dipole potential, and the slopes of the flat portions 
of E R on either side of the peak are proportional to 6, 
the derivative of the inhomogeneous magnetic field. 



H R = U + H! 



f(r) 





-2 



4 




(26) 
(27) 





" 2fi(B + bZ) 





o ■ 


+ 


-2fj,(B + bZ) 





















-gz 


_ 



where g = bfi/2 is introduced. 

Next this reduced Hamiltonian can be represented in 
the 1 5), \T) subspace as 



E- ( !lf 2j) 




H R 



4/(r) 



" 1 


0" 




' 


f " 








- gz 


l 






We express it in terms of the Pauli matrices 
H R = 2f(r)I + /3-a 



(28) 



(29) 



where I is the identity operator and (3 = (—gz, 0, 2f(r)). 
Then the eigenvalues are given in terms of the angle ui 
which is defined by 



tan lj 



-gz 



2/ 



and we choose the branch 
-gz 



-bzr 3 
~4^~ 



-bzr 3. 



smw 



VW 2 



b 2 7 



(30) 



(31) 



Physically tan u> represents the ratio of the energy of the 
dipole in the external inhomogeneous magnetic field to 
the energy due to the dipole-dipole coupling. Solving for 
the eigenvectors and eigenvalues of H 2 , we obtain 



and 



|-o) 
E R 

\+a) 



E* 



-sin||T> + cos||S) 
2/-| j 9|=2/(l-secw) 



cos % \T) 



flllr \S) 



2 1 ' 2 
2/ + |/3| = 2/(1 + sec w). 



(32) 
(33) 

(34) 
(35) 



FIG. 1: E^(zi,Z2) for an inhomogeneous field over all 
space plotted for systems in which the particles are 
on opposite sides of the origin. The system could 
be represented by a point in the near quadrant, in 
which case z\ > and z 2 < 0. If particle 1 and par- 
ticle 2 are interchanged, then the system would be 
represented by a point in the far quadrant obtained 
by a reflection across the line zi = 22- For El the 
energy is monotonically decreasing as either particle 
moves away from the origin and there is no stable 
state. All the energies are scaled by Eq = /j,Bo and 
the coordinates are scaled by ro which is defined by 
E = 2/(r ) r = (2^/Bq) 1 / 3 . 

It is interesting to determine the nature of the spin 
wavefunctions for |±a) for large and small values of sep- 
aration r when a magnetic field is present. For large 
r(z > 0) 4=> r — > 00 (sinw — > —1): 



\-a) 
E R 



1 

72 



(\T) + \S)) = \+) 1 \-) 2 , 



-gr, 

7=2 m 



\S)) 



->ll+> s 



E* 



+gr 

and for small r^r^O (sin u> 
\-a) - El 



0): 



16 r 



(36) 
(37) 
(38) 
(39) 

(40) 
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E- 




FIG. 2: E_(zi, —Zi) for an inhomogeneous field over all 
space. This is the same as the contour with 2:1 + 22 =0 
in the FIG. 1. 



E+ j 




FIG. 3: The positive eigenvalue E+(z\,Zvl) for an inho- 
mogeneous field over all space. We see stable states 
can occur in this case. 



\+a) - \T) , E\ 



(41) 



For large separations, the energy eigenvalue is dominated 
by the effect of the inhomogeneous field, whereas for 
small separations, the eigenvalue is dominated by the 
dipole-dipole interaction. As the relative position z of 
two particles goes from —00 to to +00, the state corre- 
sponding to |+a) goes from an unentangled state, namely 
|— )i |+) 2 , to an entangled triplet \T) state, and then 
back to the unentangled |— ) 1 |+) 2 - By superposition of 
states with different momenta, we could form a state that 
would oscillate in time between entangled and unentan- 
gled states [llj. 

We need to avoid energy levels crossing each other dur- 
ing the application of the adiabatic perturbation. Other- 
wise transitions between the levels may occur. From the 
reduced Hamiltonian 126|) one can easily find the energy 
levels for the other members of the triplets |2±i) to be 
— 2/(r) ± 2//(£>o + bZ). Therefore, under the assumption 
that the two particles are on the opposite sides of the ori- 
gin {Z ps 0), we require the following equalities to avoid 
crossing of energy levels: 



2/(r) 



2iiB > E* > E R > -2/(r) - 2 M B , (42) 



which are obeyed provided the energy contribution from 
the constant magnetic field Bq is significantly greater 
than the contributions from either the inhomogeneous 
field or the dipole potential (fiB > 3/(r) + g\z\/2) [l^ . 

Independently we also need a condition in order for the 
adiabatic evolution to proceed (see e.g. Q). Namely we 
need to turn on the magnetic field (inhomogeneous part) 
slowly, and the condition for the time period T for this 
switching is 



Kgzf 



(4/)VCff*) 2 + (2/) s 



(43) 



Et(:i f -zi) 

2Qr 




FIG. 4: E + (zi, —Zi) for an inhomogeneous held over all 
space. The rapid increase is due to the divergence of 
the dipole-dipole potential. 



32^ 3 y/b 2 r s + 16^ 2 



(44) 



Therefore the greater the separation, the more slowly the 
inhomogeneous field needs to be applied. 



A. Use of Born-Oppenheimer Approximation to 
the Tunneling effect 

So far we have neglected the kinetic term under the as- 
sumptions that the characteristic frequencies of the spin 
precession are much greater than those of the transla- 
tional motion of the two particles. These assumptions 
are strictly true, for instance, for the NMR case where 
particles are part of the molecules and hence they are 
always bound. 

In contrast, we want to consider a situation where the 
particles are essentially free. For the eigenstate \+a), a 
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graph of the corresponding eigenvalue E+(z) is a pos- 
itive double hump function which describes two bound 
particles, one on either side of the origin. The potential 
permits tunneling across the barrier. In this tunnelling 
process the two particles would be exchanged. First we 
note that the Schrodinger equation for the center of mass 
motion as well as the relative motion is invariant under 
the transformation z\ — > z 2 and z 2 — > z\ which corre- 
sponds to the tunneling transition. In (zt, z 2 ) space, this 
transformation corresponds to a reflection across the line 
z\ = z 2 . The corresponding states are degenerate in en- 
ergy. The Hamiltonian is also invariant under the parity 
transformation z\ — ► —z%, z 2 — > — z%. 

In order to analyze this process, we can utilize a Born- 
Oppenheimer approximation and we can separate two 
degrees of freedom by the use of an average potential in 
the Schrodinger equation for the relative motion of the 
particles. Thus the problem is reduced to a one body 
problem with the potential for the relative coordinates. 
The appropriate potential is in fact the energy eigenvalue 
E+(z) which was obtained in the previous section: 



Therefore this approximation yields: 



H BO \E) = E\E) 
P 2 



H 



BO 



+ E*(z). 



(45) 

(46) 
(47) 



As is usual for the singular potential oc r~ n (n > 2) the 
barrier becomes infinitely high at the origin. To estimate 
the tunneling probability we introduce a cutoff for small 
r defined by r c . Figure 5 shows a typical shape of the 
potential barrier after this regularization. 



(Ea) is the eigenvalue of \<ps) (10a))- Now we consider a 
state in which the particle is located on either the right 
side of the barrier or the left side and express these states 
by \4>r) and \4>l) respectively. Notice that they are not 
eigenstates of the Hamiltonian (|47|l in general. Obviously 
we have the following relation: 



1 



v/2 



(\<Pr) + 



\4>a) = -^(|^>-|<M) 



(48) 
(49) 



up to a phase factor. Therefore if we choose the initial 
state at t = as \4>r), then after simple calculations we 
get the state at time t as follows. 



n 



•isin(^)|<fo) 



(50) 



where A = (Ea — Es)/2 and we omit an overall time 
dependent phase factor which is irrelevant right now. 
Therefore we observe that the state is oscillating between 
a configuration in which the particle is on the right side 
and a configuration in which the particle is on the left 
side. In terms of the original variables, the positions of 
two particles are being exchanged in time. 

To know the characteristic frequency A/h we need to 
know the eigenvalues and the states of the Hamiltonian 
using the relation A = — (4>r\Hbo\<Pl)- This requires 
knowledge about the solution to the Schrodinger equa- 
tion 146|) . Instead we can use the WKB approximation 
to examine the possibility of this effect. 

Using the WKB method we calculate the probability 
w for a particle at rest around one minima r m of the 
potential to tunnel across the potential barrier to the 
other minimum: 




FIG. 5: Potential energy barrier E+(z 
the regularization. 



Z2) after 



Since the Hamiltonian H47[) commutes with the parity 
operator, the eigenstates of it are also eigenstates of the 
parity, namely they are either symmetric or antisymmet- 
ric. Let us denote these states by \<j>s) and \4>a) respec- 
tively. Then in general we know E$ < Ea, where Es 



)(E = 0) ~ exp 



m/i 2 . 3 
h 2 VJV ^ ' 



(51) 



A cutoff distance of 10 _15 m was chosen for r c , approx- 
imately the Compton wavelength of the neutron. This 
value for r c and the value of r m from Appendix B yields 
an estimate for w ~ e -°- 94 ^ 0.39. It suggests a possibil- 
ity of the tunneling effect, namely, the exchanging of the 
two particles through the barrier. We summarize these 
calculations in Appendix B. 

Two remarks are in order: when the potential bar- 
rier becomes infinite the ground state will be degenerate 
(Es = Ea <=>• A = 0). This means there are no oscilla- 
tions at all. And hence there exist the states \<t>R) and 
\4>l) separately. Although these states are not eigenstates 
of the parity operator, because of the degeneracy they are 
allowed states. This is an example of the spontaneous 
symmetry breaking. To examine the possibility of the 
exchanging effect more precisely we also need to include 
the hyperfine interaction term in the original potential 
P)l. which becomes important at short distances. 
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IV. CASE 2: MAGNETIC FIELD PRESENT ON 
THE RIGHT SIDE OF THE ORIGIN ONLY 

We assume the magnetic field B(z) is in the z-direction 
and that it is non-zero only in some region on the right 
side of the origin (z\ ,z 2 >0). Outside of this region, for 
example on the left, the magnetic field vanishes so, for 
example, B(zi)iPl(zi) = 0, B(z 2 )iPl{z 2 ) = 0, etc. This 
magnetic field breaks the translational symmetry of the 
field present in Case 1. Using these properties we can 
show that 



E+ \S) = E+ \T) = 

£_|S) = -B T \T) 
E_|T) = -B T \S) 



where we define 



B T (z u z 2 ) = B(zi) + B(z 2 ) 



(52) 

(53) 
(54) 



(55) 



Note that for our model, either particle 1 or 2 will be on 
the left side of the origin so either B{z\) or B(z 2 ) will 
vanish. Thus for the special case of a constant field Bo on 
the right side, we have Bt = Bq. With our assumptions 
we can express Bt{z\ 1 z 2 ) in the following way: 

B T ( Zl ,z 2 ) = B[z 1 )Q{z 1 - z 2 ) + B{z 2 )9(z 2 - Zl ) (56) 

In the impulsive approximation, these results yield a 
representation of H] using the basis states {\T) , \S)} 



Hi = ixB T 



1 

1 



(57) 



We notice that the Bt field causes transitions between 
the singlet and triplet states like the inhomogeneous field 
in all space. This suggests that even the simplest case (a 
constant field Bq on the right side) will do this. This 
transition might be of interest, in particular, in applica- 
tions of quantum computation. Similar discussion were 
done by several authors (see Ref. Q and references 
therein). However, the transition between the singlet and 
the triplet was not discussed in there. 

The remaining components of the triplet also trans- 
form among themselves, which follows from: 

£+|T ±1 ) = ±Bt|T±i) (58) 
E_|T ±1 ) = 0. (59) 

The interaction Hamiltonian in the basis , } 



Hr 



fiB T 



1 

-1 



(60) 



In terms of the basis {|T_i) , |Ji) , \T) , \S)} , the total 
reduced Hamiltonian H R can therefore be written as 



H R = f(r 
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liBi 
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0-100 
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10 



(61) 



After applying the magnetic field perturbation, the 
|T_i) and |Ti) components still do not mix with any other 
components. The total reduced Hamiltonian is diagonal 
in the subspace spanned by |T_i) and \T±) , and the cor- 
responding energy eigenvalues are E R 1 = —2f+fj,BT and 
E R = ~2f-(iB T . 

We now consider in detail the subspace spanned by 
|T) , \S) which is not diagonal. In the same manner as in 
the previous section we express the total reduced Hamil- 



tonian H R in terms of the Pauli matrices as: 



Hi 



2/(r)J 



(62) 



where the vector a. = (fiBx, 0, 2/). Define the angle 9 by 

a x fj,B T (zi,z 2 ) B T r 3 



tan 8 

a z 2f 

and the quadrant is specified by 

fiB T 



2fi 



sin# 



^bJTW' 



(63) 



(64) 



Solving for the eigenvectors and eigenvalues of H R , we 
obtain 



\-a) = - s in-|T>+cos~|S>, 
E R = 2/- \a\ = 2/(1 -sec 0) 



and 



l+«> 
E R 



cos-|T)+sin-|S), 
2f+\a\ = 2/(1 + sec ( 



(65) 
(66) 

(67) 
(68) 



To avoid the crossing of energy levels, we require E R > 
E R 1 which implies that 

2(mB t < 3/(r). (69) 

And also from the adiabatic theorem we get 



h(fiB 7 



(4/) V(^t) 2 + (2/) 5 



(70) 



Using the above condition l|69|l we estimate the lower 
limit of the time T m i n as 



9 % 

so 7M 



0.11^/" 



(71) 



The restriction of Eq. I|69|) severely limits possible 
modulation frequencies of the field being applied adia- 
batically. We briefly mention a procedure, similar to that 
used in Case 1, which results in a much larger bandwidth. 
We first apply the large constant B field impulsively over 
all space [l4|. Then we apply the inhomogeneous field 
adiabatically over the right side only. With this proce- 
dure, the restriction on the energy levels to avoid cross- 
ings is similar to Eq. lj4*2)l . but the levels |T_i) and 
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are also shifted by the inhomogeneous field ±/j,Bt\b =o- 
The quantity hBt\b =o i s always positive since the inho- 
mogeneous field is on the right side, so the requirement 
for no level crossing is less stringent, and is met simply 
if fiBo > 4/. The corresponding requirement for the va- 
lidity of the adiabatic approximation is the same as Eq. 
(I44|l . with r replaced by the positive value of the pair 
{zi,z 2 }. 



A. Homogeneous field on the right side only 

We can consider several different magnetic field 
strengths Bt{zi, z 2 ). For the simplest case, in which Bt 
is a constant on the right side, we find that E^(zi,z 2 ) 
is a negative potential (Figures 6-7) that does not con- 
fine the particles so no stable states are expected, while 
E+(zi,z 2 ) is a positive potential that also does not con- 
fine particles (Figures 8-9). 




E- (s L ,-si) 




-0 ,8 



-1 



FIG. 7: E-(zl, — zl) for a constant magnetic field This 
curve is obtained by plotting the function E-(zl,z1) 
for the contour zl + z2 = . 




FIG. 6: E_(zi,Z2) for a constant magnetic field on the 
right side, plotted for particles on opposite sides of 
the origin. 



FIG. 8: E+(zl, z2) for a constant magnetic field on the 
right, plotted only for particles on opposite sides of 
the origin. 



B. Inhomogeneous field on the right side only 

Next we shall consider an inhomogeneous field (Bt = 
Bo + bZ + br/2) on the right side. The graphs for 
E^(z\,z 2 ) (FIG. 10-11) are plotted with the same pa- 
rameters as FIG. 1-4 for the comparison. Similar results 
are obtained in this case. 

The positive energy E^(zi,z 2 ) has a double hump 
structure (FIG. 12-13) with two potential wells for quasi- 
bound states separated by an energy peak. We note that 
unlike in case 1, there is no symmetry with respect to 
reflection across the line z\ + z 2 = so the forces act- 
ing on the particles are not equal. The force due to the 
magnetic field acts only on the particle on the right side. 




FIG. 9: E+(zl, — zl) plotted for a constant magnetic 
field on the right side of the origin. 



Here we define the effective forces through partial differ- 
entiations of the effective potentials E± with respect to 
z\ and z 2 - 




E- (s L , z 2 y 




15\ 

E+(3 lf Sl ) 1 




FIG. 12: E+{zl,z2) for an inhomogeneous field on the 
right side, plotted for particles on opposite sides of 
the origin. 



FIG. 10: E-(zl, z2) for an inhomogeneous field on the 
right side, plotted for particles on opposite sides of 
the origin. All parameters for the FIG. 10—13 are 
chosen same as for FIG. 1—4. 




FIG. 11: E-(zl, —zl) for an inhomogeneous field on the 
right side only. 



C. Summary 

In summary, with the adiabatic evolution of spins when 
we apply the magnetic field on the right side, we find that 
the initial singlet ground state singlet \S) has evolved into 
the eigenstate given by a linear combination of singlet 
and triplet states: 

|-o) = -sin^|T) + cos^|S). (72) 

In this adiabatic evolution, the reduced energy goes from 
for \S) to 2/(1 — sec9) for |— a). Since \—a) is a spin 



eigenstate of the total Hamiltonian, there is an over- 
all time dependent phase factor, e ~ lE - t / h for this state, 
which we have omitted. 

The eigenstates |±a) have the property that, for any 
value of 9, the spins of the particles in the z-direction are 
always correlated: 



(S Zl +S Z2 )\±a)=0. (73) 



This result, which follows since [Ht, S Zl + S Z2 ] = 0, indi- 
cates that the spins of the two particles remain correlated 
if an inhomogeneous magnetic field is applied only to the 
particle on the right side of the origin. 
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V. COMPUTATION OF THE DENSITY 
MATRIX FOR AN ADIABATIC CHANGE IN 
THE MAGNETIC FIELD ON THE RIGHT SIDE 
ONLY 

In this section we shall study the case 2 further in de- 
tail. One convenient way to consider the effect of the 
magnetic field (which is present only on the right side) 
on the state of the system on the left side (where the 
magnetic field vanishes) is to predict the results of a mea- 
surement of an observable A which is defined only on the 
left side. The predictions can be done by means of the 
reduced density matrix p(L, —a) on the left side L for the 
state | -a) El: 



p(L, -a) = Trright side | -a) (-a| 



(74) 



In order to compute this reduced density matrix for |— a), 
we first compute the reduced density matrices for the 
basis states \S) and \T). For an observable A which is 
defined only on the left side we have Atpn(z 2 ) = 0, etc. 
The result of a measurement of A for the state |T) is 
defined as 

Trp(L, T)A = [T\ A \T) (75) 
and can be computed using the expression for |T): 

(T\A\T) = l[ 1 (+L\A\+L) 1 + 1 (-L\A\-L) 1 (76) 
+ 2 (-L\A\ - L} 2 + 2(+L\A\ + L) 2 ] (77) 
where we have defined 



\±L) 
\±L) 



l±>i^z(zi) 
l±) 2 ^(*a)- 



(78) 
(79) 



The corresponding reduced density matrix on the left side 
is 



P(L,T) 



1 



-L\ 



-U 



(80) 



+ \ + L) 22 (+L\+ \-L) 22 (~L\}. (81) 



By a similar calculation we find that the reduced den- 
sity matrix for \S) has the same value as for \T): 



p(L, S) = Tr rig htside \S) (S\ 
= P(L,T). 



We also compute 

p(L, TS) = Tr r i g htside 



\T)(S\ 



(82) 
(83) 



(84) 



= -[|+L) u (+L|-|-L) n <-L| (85) 
+ \+L) 22 (+L\-\-L) 22 (-L\] (86) 



and find that 



Using the preceding results, we can write an expression 
for p(L,—a) using the definition of the state |— a): 



p(L,-a) = sin 2 -p(L,T) + cos 2 -p(L,S) 

-2sin^cos^p(L, TS) 
= p(L, S) — p(L, TS) sin 6. 



(88) 
(89) 



This last result can be written in matrix form in the 
basis {\+L) 1 ,\~L) 1 ,\+L) 2 ,\-L) 2 \ : 



p(L,-a) = - 



1 



- sin ( 








- sin ( 








- sin ( 








1 + sin ( 



(90) 

One can contemplate making standard quantum me- 
chanical measurements in which a superposition collapses 
to a eigenstate0|- Alternatively we consider the use of 
a protective Stern-Gerlach measurement, in which there 
is an adiabatic interaction between the pointer and the 

system[lH[l3. 

We want to examine the possibility of 
the protective measurement in our model under the as- 
sumption we will work in a decoherence free subspace. 
In previous discussions of protective measurements, only 
the case of a small perturbation was considered in or- 
der to minimize the change in the wavefunction of the 
state to be measured. In our system, such restrictions 
are not necessary. As long as the interaction is adia- 
batic, the state will evolve continuously as an eigenstate 
of the instantaneous Hamiltonian, without any transition 
provided there is no degeneracy in the energy. If the mag- 
netic field is applied adiabatically, the field can become 
large, resulting in a large change in the wavefunction. 
Because the changes are adiabatic, they are reversible as 
the magnetic field is reduced. Since a protective mea- 
surement of the spin for our two particle system does not 
change the state of the system, the protective measure- 
ment does not end the entanglement. 

In the spin-spin model, the states |— a) and |+a) into 
which | S) and |T) evolve when a magnetic field is ap- 
plied on the right side are non-degenerate so a protective 
measurement of the reduced density matrix should be 
possible. Indeed all four basis states|— a), |+a), |Ti), and 
|T_i) are non-degenerate provided / ^ and B ^ so 
a protective measurement should be possible of the den- 
sity matrix[6j . As a consequence the result of a protective 
measurement of an observable A that acts only on the left 
side can be written as the trace over the reduced density 
matrix: 



< A > {L _ a) =Tr(p(L,-a)A). 



(91) 



p{L,TS) = p{L,ST). 



(87) 



Note that here < A > represents the result of a single 
protective measurement on a single system. It does not 
have the usual meaning of the expectation value of A, 
which is based on the measurement of the observable A 
for an ensemble of identically prepared systems. Since 
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the elements of the density matrix depend on 9, it is clear 
that changes in 9, resulting from either changes in B or 
(z\ — Z2), will affect the measured values of obscrvables 
A. In other words, in this model with a finite separation, 
by changing the value of the magnetic field on the right 
side, it is possible to detect the effect on the left side 
nonlocally. Note that nowhere do we use the specific 
form of the potential. Indeed, we could use the formalism 
presented with any function of the separation. 

To illustrate the application of a protective measure- 
ment, we could measure the total spin in the z- direction 
for both particles using local measurements on the left 
side only. We assume we use a Stern-Gerlach analyzer 
with an inhomogeneous field in the z-direction and we 
observe the deflection of the wave packet corresponding 
to a particle. We assume we have calibrated our appa- 
ratus so that we can determine the spin of a particle on 
the left side by observation of the deflection in the z- di- 
rection at a certain point in the experiment. Since the 
particles are indistinguishable, a single protective mea- 
surement of the spin of a particle on the left will yield a 
deflection corresponding to the result 

(S Z1 + S Z2 ) left Slde = tr{p[S Zl + S Z2 ]) = ~ sin0 (92) 

where the value of 9 is determined by the magnetic field 
on the right side and the distance between the particles. 
The corresponding total spin in the z-direction measured 
on the right side is isin#. Since this protective mea- 
surement has not collapsed the wave function, we could 
do additional measurements on the same state to de- 
termine the rest of the density matrix P(L-a) in the 
\T) ,\S) subspace. If we used Stern-Gerlach analyzers 
with fields in the x and y directions for two additional 
experiments, respectively, we would find: 

(S Xl + S X2 ) left side = (S Vl + Sy 2 ) left side = 0. (93) 

These components of the spin vanish because we chose 
the magnetic field on the right side to be in the z- 
direction. 

It is interesting to contrast the protective measurement 
of S Zl +S Z2 with the ordinary measurements. In an stan- 
dard quantum mechanical measurement the wavefunc- 
tion collapses and an eigenvalue of the operator is mea- 
sured. To determine the expectation value of the opera- 
tor, one needs to make a statistically significant number 
of measurements on an ensemble of identically prepared 
systems. In our model, since the particles are indistin- 
guishable, for each standard measurement on a different 
identically prepared system, a deflection in the Stern- 
Gerlach apparatus will be measured that corresponds to 
spin of +i or — i. If enough measurements are done, 
it will be determined that the probability of measuring 
+ ^ is |(1 — sin 9) and the probability of measuring — \ 
is 5(1 + sin 9). Using these probabilities, the standard 
measurements will determine that the expectation value 
of the operator S Zl + S Z2 on the left side is — ^sin^, in 
agreement with the single protective measurement. 



This comparison of the two methods illustrates some 
of the advantages of protective measurements and some 
of the disadvantages of standard measurements. 

VI. CONCLUSION AND DISCUSSIONS 

A system was considered in which two neutral spin 
1/2 particles interact through a diple-dipole potential. 
The potential leads to singlet and triplet entangled states 
which are modified when we apply a magnetic field over 
all space (Case 1), or over just the region to the right of 
the origin (Case 2). 

For the case 1 we showed that the singlet state is es- 
sentially unstable and it tends to separate. However, the 
spin zero component of the triplet will have a double 
hump shape energy which implies the particles will tend 
to stay near the minima with high probabilities and form 
a bound state. Also we showed that it is possible to have 
a tunneling effect in this case. Such an effect might have 
applications in quantum computation. 

For the case 2 we observed that although the external 
field is applied only in one part of the system, the other 
part will be affected by it due to the entanglement of the 
system. This manifests the nonlocality of the entangle- 
ment in quantum mechanics. It appears protective mea- 
surements can be used to determine the density matrix 
without ending the entanglement. As discussed above for 
the case of finite separation, one limitation in the signal- 
ing between the two regions of space in this model lies 
in the requirement that we must have an adiabatic per- 
turbation or transitions will occur between states. When 
the magnetic field is turned on, it must be done slowly 
enough so that no transitions are induced between the 
initial state and other states. Although the specific re- 
strictions depend on the manner in which the perturba- 
tions are applied, they require that the potential due to 
the dipolc-dipole interaction does not vanish. 
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APPENDIX A: MODEL WITH CONFINING 
POTENTIALS 

In this appendix we will examine our model with ad- 
ditional confining potentials (for instance, created by the 
optical laser). For simplicity we approximate this po- 
tential by a harmonic potential around certain point 
Xq with characteristic frequency 0, namely, V(x) — 
^mfl 2 (x — Xo) 2 . We choose the center of the potential 
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in such a way that two particles will be separated by a 
distance z . Under these conditions the ground state en- 
ergy of the particles is §fi (we use K = 1 in appendices.) 
and the corresponding states are 

Mxi) = (i)i e -i«*?-^+(»«-*)") (Al) 

7T 

^( Xi ) = (i)f e -|S(^+!/?+(-i+^) 2 ) ; (A2) 

for j = 1,2 and the parameter £ = \/ mil is introduced. 
To have a stable ground state in our model we require 
the energy levels are well separated, or the parameter £ 
is large. (Similarly we can assign these wave packets for 
the free particles case too.) After we obtain the spatial 
wave functions for the singlet and the triplets [18j , we can 
estimate the effect of the hyperfine interaction term in the 
dipole-dipole potential ©■ We evaluate the expectation 
values of the delta function as follows. 



(5(x)) s = 4(|-)§(1 
(S(x)) T = 0. " 



(A3) 
(A4) 



Therefore, under our assumptions the contribution from 
the hyperfine interaction is negligible in our model be- 
cause of the exponential factor in (|A3|> . However, in con- 
trast, this term plays an important role in a situation 
where the distance becomes very small and we need to 
add this contribution too. This regime is discussed in 
the literature, for instance, see Ref. and references 
therein. 

Next we estimate the change in the positions of the 
particles and the kinetic energy of them. From simple 
calculations it is easy to see that the expectation values 
of the relative position and the center of mass are zero. 
Similarly we obtain for the expectation value for the cen- 
ter of mass kinetic energy 



p2 

2m 



p2 

2m 



2m 



S \ """I T 

and for the relative kinetic energy we have 



P_ 

m, s 
m 



4m 
4m 



(1 



(1 



t4 



2 2 £ z o 4- 1 ' 



t4 



2 2 £ Z Q — 1 ' 



(A5) 

(A6) 
(A7) 



Again, the second term inside the parenthesis is negli- 
gible by our assumptions. And hence we see that the 
corrections due to the motion of the particles are small. 



APPENDIX B: CALCULATION OF THE 
TUNNELING PROBABILITY (tSTll 

To estimate the tunneling probability we simplify our 
model as follows. First we calculate the minima r m of 
the potential jjSjl , this is easily found 



= 2VT5f(r rn )^r m = (^-)-s. (Bl) 



Then we evaluate the the tunneling probability between 
two minima through the dipole-dipole interaction barrier 
assuming that the potential is mainly dominated by this 
interaction in small distance. We use the WKB method 
[l3| to get the probability w(E) for the particle with the 
kinetic energy E = k 2 /m (the mass is m/2), that is, 



w(E) ~ exp[-2 



y/m{V{z) - E) dz 



(B2) 



In our case the simplified potential is V(z) = 4/i 2 /|z| 3 . 
As we mentioned earlier we introduce the cutoff at r c 
to get rid of the divergent integral. Under these as- 
sumptions and defining the momentum k m and k c corre- 
sponding to the minima and the cutoff respectively, i.e. 
km/ m = V{ r m) and k 2 /m = V(r c ), the integral inside 
the exponent of i|B2|l is evaluated as below. 



2^kJ^r m -kr m (^B L ,( — ) + 3VkJ^r. 



(B3) 

where B a ^(x) is the incomplete beta function defined by 
B a ,b(x) = Jq t a ~ 1 (l — t) b ~ 1 dt and we neglect the positive 
orders of r c (in r c /r m <C 1 limit). Therefore for the rest 
particle (k = 0) at the minima, a rough estimate for the 
tunneling probability is 

w{E = 0) ~ exp[-4(3fc c r c - 2k m r m )] (B4) 

which is given in the text after substituting k m and k c . 
Another case is for the particle having the minimum po- 
tential energy yk — km)'. 



w(E = V{r m )) ~ exp[-4(3A: c r c - 5(1 \)k m r m )] (B5) 

o 2 



where £?(§, |) — 2.24 is the usual beta function. And 
hence we see that this is the same order as (|B4[1 . 
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